Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1659 1651 (N30W72 L= 53)   976     77        1 (1616/1633-1701/1671) -21.8 W66-W78 N27-N32
      2 1383 1542 (N31W36 L= 18)   458     75        1 (1348/1518-1430/1560)  -8.1 W33-W40 N29-N32
      3  668 1330 (N30E19 L=323)  1930    523        3 ( 458/1247- 891/1439) -25.9 E37-E04 N22-N40
      4 1116 1132 (N11W07 L=349)   219     49        1 (1092/1106-1138/1163)  33.7 W05-W08 N10-N13
      5  440 1036 (N15E38 L=304)   334    138        0 ( 408/ 985- 500/1069)  29.7 E40-E33 N12-N16
      6  869 1016 (N08E09 L=332)   229     62        1 ( 852/ 978- 890/1048)  45.0 E10-E07 N06-N10
      7  285  931 (N10E53 L=289)  1526    336        3 ( 200/ 791- 369/1040) -50.9 E61-E46 N01-N17
      8 1526  941 (S07W30 L= 11)   227     85        1 (1513/ 911-1551/ 987)  63.4 W28-W31 S08-S04
      9 1171  790 (S10W05 L=347)   665    126        2 (1151/ 728-1194/ 865)  69.4 W02-W06 S13-S05
     10 1393  731 (S17W18 L=360)   446    176        0 (1356/ 701-1454/ 769) -23.2 W15-W22 S19-S14
     11  653  599 (S14E31 L=311)   809    287        0 ( 599/ 541- 725/ 658) -26.6 E33-E25 S17-S09
     12  872  392 (S31E23 L=319)  1117    243        2 ( 787/ 330- 956/ 488)  23.6 E30-E14 S35-S25
  

    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