Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  737 1790 (N44E45 L=209)   619     79        1 ( 708/1755- 768/1824) -24.4 E50-E39 N41-N46
      2 1663 1506 (N42W44 L=298)   639    131        2 (1602/1465-1719/1536)  -4.4 W36-W50 N39-N43
      3  418 1484 (N17E52 L=203)   590    133        1 ( 366/1450- 472/1506)  26.6 E55-E47 N13-N18
      4  931 1352 (N17E11 L=243)   196     62        2 ( 905/1329- 960/1371) -14.0 E13-E09 N16-N18
      5  879 1267 (N11E13 L=242)  1746    299       12 ( 782/1193- 994/1350)   4.1 E19-E05 N06-N16
      6  603 1281 (N07E30 L=224)   529    143        0 ( 586/1208- 621/1345) -53.1 E33-E27 N03-N11
      7  791 1244 (N08E18 L=237)   311     79        1 ( 762/1206- 815/1277) -36.9 E19-E15 N06-N09
      8 1068 1129 (N06W01 L=255)   255    114        2 (1031/1093-1120/1146)   0.0 E01-W04 N03-N07
      9  843 1108 (N01E12 L=243)   199    116        0 ( 809/1060- 881/1142)  59.0 E13-E10 S02-N03
     10  371 1082 (S08E44 L=211)   406     93        1 ( 324/1063- 422/1103)  20.6 E47-E39 S09-S06
     11  842  890 (S12E09 L=246)  1207    379        1 ( 764/ 746- 940/1050)  74.1 E11-E05 S21-N00
     12  414  777 (S26E38 L=217)   703    183        1 ( 335/ 744- 471/ 835)  30.3 E45-E32 S28-S21
     13 1305  714 (S14W23 L=277)   385    138        0 (1241/ 674-1358/ 748) -12.5 W18-W27 S15-S13
     14 1061  662 (S22W09 L=263)   346    106        1 (1014/ 640-1117/ 685)  38.7 W05-W12 S23-S19
     15 1307  581 (S22W27 L=281)   320    107        0 (1256/ 564-1358/ 604)  26.6 W23-W29 S24-S20
     16 1138  502 (S31W18 L=273)   319    122        1 (1085/ 479-1190/ 546)  51.3 W14-W21 S33-S26
     17  384  429 (S53E53 L=201)  1247     93        0 ( 352/ 398- 417/ 467)  28.3 E61-E45 S56-S49
     18  469  379 (S56E42 L=212)   845     72        1 ( 439/ 342- 498/ 421)  81.9 E47-E38 S58-S51
  

    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