Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1502 1613 (N24W49 L= 84)  3428    581        4 (1264/1479-1740/1681) -33.7 W31-W62 N10-N31
      2 1714 1614 (N18W70 L=105)  3185    216        4 (1615/1547-1801/1677) -53.1 W62-W82 N12-N24
      3  670 1573 (N45E10 L= 24)   268     38        1 ( 643/1556- 697/1589) -26.6 E12-E07 N43-N46
      4 1389 1504 (N21W35 L= 70)   263     87        0 (1344/1491-1436/1515) -31.0 W32-W37 N19-N22
      5 1728 1213 (S03W49 L= 83)   310     53        0 (1711/1183-1747/1239)  71.6 W47-W50 S04-S01
      6 1730 1117 (S09W46 L= 81)  1504    359        3 (1639/1046-1792/1183) -36.0 W39-W51 S12-S03
      7  442 1134 (N22E34 L=  1)   749    189        1 ( 353/1075- 505/1179) -39.3 E40-E29 N17-N26
      8  915  865 (S03E10 L= 25)  1849    489        5 ( 774/ 685-1056/1002) -62.2 E15-E05 S15-N05
      9 1347  908 (S11W16 L= 50)   835    258        3 (1267/ 794-1390/ 992)  38.7 W08-W19 S15-S05
     10 1221  897 (S08W08 L= 43)   434    130        3 (1191/ 833-1242/ 961)  54.5 W05-W10 S11-S04
     11 1272  928 (S08W12 L= 46)   236     48        1 (1254/ 903-1297/ 954)  45.0 W10-W13 S09-S06
     12  382  904 (N10E42 L=352)   607    146        1 ( 324/ 873- 451/ 946) -29.1 E46-E37 N07-N13
     13 1475  676 (S28W19 L= 54)  1158    270        4 (1383/ 635-1579/ 724) -23.2 W13-W27 S29-S23
     14  587  677 (S07E34 L=  1)   255     34        1 ( 564/ 651- 608/ 697) -33.7 E35-E32 S07-S05
     15 1156  410 (S37E08 L= 26)  3654    686        6 (1005/ 196-1258/ 695)  45.9 E30-W02 S49-S18
     16  572  600 (S11E38 L=357)   727    143        3 ( 512/ 562- 647/ 634) -26.6 E41-E32 S13-S07
     17  538  419 (S21E49 L=346)  2868    478        6 ( 389/ 316- 656/ 562)  13.1 E77-E34 S24-S14
     18 1068  338 (S39E17 L= 17)   247     76        0 (1061/ 294-1079/ 380)  50.2 E19-E14 S42-S36
  

    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