Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1342 1744 (N37W49 L=290)  2271    375        3 (1191/1716-1545/1771) -24.0 W36-W65 N29-N42
      2 1596 1622 (N22W59 L=300)  1595    269        2 (1472/1563-1696/1676) -38.2 W51-W67 N16-N28
      3 1563 1543 (N19W51 L=291)  1975    346        0 (1442/1428-1703/1652) -67.2 W47-W57 N08-N28
      4 1060 1551 (N33W17 L=258)   269     83        0 (1019/1539-1106/1562) -18.4 W14-W20 N31-N34
      5  640 1366 (N33E16 L=225)  1129    466        2 ( 441/1325- 841/1400) -16.7 E32-E02 N27-N36
      6 1083 1169 (N11W07 L=248)   495    157        1 (1037/1098-1124/1236) -84.3 W05-W08 N05-N15
      7  790 1200 (N19E09 L=231)   205     83        1 ( 749/1186- 836/1214) -18.4 E12-E06 N17-N20
      8 1801 1167 (S08W55 L=296)   486    102        1 (1769/1112-1830/1203)  33.7 W50-W58 S10-S06
      9 1624 1027 (S11W36 L=277)   302     72        0 (1595/ 991-1645/1064) -78.7 W34-W37 S13-S08
     10 1734 1036 (S13W45 L=286)   285     81        1 (1690/1016-1760/1064)   8.1 W41-W48 S14-S12
     11 1863  907 (S25W58 L=298)   496     83        0 (1843/ 866-1885/ 947)  45.0 W54-W60 S27-S23
     12 1726  712 (S33W41 L=282)   799    135        0 (1689/ 672-1771/ 780) -51.3 W37-W47 S37-S27
     13 1565  660 (S31W26 L=267)   744    217        3 (1497/ 598-1640/ 745)  17.1 W20-W33 S34-S28
     14  734  657 (S10E26 L=215)   223    117        0 ( 722/ 600- 743/ 713)  66.8 E27-E24 S13-S06
     15  532  642 (S07E39 L=201)   256    101        0 ( 480/ 626- 582/ 654)  -8.1 E43-E36 S08-S05
     16  991  552 (S23E14 L=227)   800    151        1 ( 955/ 490-1012/ 632)  52.1 E17-E10 S26-S17
     17 1426  454 (S42W12 L=253)  1123    241        0 (1322/ 382-1491/ 532)  25.0 W02-W18 S44-S37
     18 1031  199 (S50E31 L=210)  1769    264        3 ( 904/ 156-1141/ 253)  -2.2 E45-E18 S52-S47
  

    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