Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1538 1743 (N31W65 L=103)   615     66        1 (1515/1713-1557/1776)  29.7 W61-W71 N28-N32
      2 1653 1609 (N20W63 L=101)  1740    178        0 (1585/1549-1722/1663) -68.2 W58-W68 N14-N24
      3 1467 1583 (N24W45 L= 83)  3137    661        7 (1213/1400-1698/1656) -41.5 W27-W58 N07-N31
      4 1335 1480 (N21W31 L= 69)   257    111        1 (1283/1463-1389/1491) -26.6 W27-W33 N19-N22
      5 1840 1129 (S11W57 L= 96)  1100    175        2 (1792/1063-1883/1182) -69.4 W53-W63 S14-S06
      6 1720 1089 (S10W45 L= 83)   944    202        4 (1670/1025-1751/1159)  60.3 W39-W48 S13-S05
      7  409 1121 (N22E37 L=  1)   571    136        1 ( 351/1069- 451/1155) -45.0 E41-E34 N18-N25
      8 1385  952 (S09W19 L= 57)   273    103        0 (1366/ 901-1401/1000)  80.5 W17-W19 S11-S05
      9  857  858 (S02E13 L= 25)  1379    378        4 ( 744/ 692- 938/ 972) -80.5 E17-E10 S13-N05
     10 1298  906 (S10W13 L= 51)   601     96        3 (1245/ 873-1336/ 948)   0.0 W09-W15 S11-S07
     11 1221  906 (S08W08 L= 46)   216     36        2 (1202/ 880-1240/ 933)  56.3 W07-W09 S09-S06
     12  337  893 (N11E47 L=351)   553    117        1 ( 288/ 862- 392/ 931) -36.9 E50-E42 N07-N13
     13 1166  855 (S09W04 L= 42)   232     87        1 (1132/ 824-1189/ 895)  36.9 W01-W05 S10-S07
     14 1422  654 (S28W16 L= 54)  1121    275        2 (1320/ 617-1525/ 709) -23.2 W09-W23 S30-S21
     15 1436  681 (S26W17 L= 55)   244     73        1 (1398/ 672-1479/ 691) -18.4 W14-W20 S27-S25
     16 1117  400 (S36E11 L= 27)  3707    655       10 ( 959/ 189-1214/ 673)  42.3 E35-W00 S49-S18
     17  560  588 (S11E39 L=359)   482    127        4 ( 526/ 556- 612/ 626) -59.0 E42-E35 S13-S08
     18  439  452 (S17E56 L=342)  1541    257        5 ( 409/ 346- 481/ 565)  39.5 E66-E47 S23-S09
     19  574  439 (S20E45 L=353)  1038    173        4 ( 524/ 381- 610/ 535)  28.3 E51-E38 S23-S14
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y