Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1249 1670 (N39W15 L=165)   727    199        0 (1169/1640-1355/1692)   3.8 W08-W23 N36-N40
      2  374 1619 (N34E54 L= 96)  2090    325        3 ( 326/1566- 435/1677) -21.8 E62-E47 N29-N39
      3  214 1365 (N17E61 L= 88)  1700    388        3 ( 133/1244- 274/1442)  50.2 E71-E54 N09-N21
      4  281 1320 (N13E53 L= 96)   612     95        4 ( 255/1279- 305/1360)  90.0 E55-E51 N10-N15
      5  750 1259 (N09E17 L=132)   344    107        0 ( 719/1212- 784/1314) -56.3 E19-E15 N06-N12
      6  657 1244 (N08E23 L=126)   287    105        0 ( 601/1228- 708/1256)  -8.1 E26-E19 N07-N09
      7  632 1093 (S01E24 L=125)   365    137        5 ( 585/1056- 674/1126)   0.0 E27-E21 S02-N01
      8  310 1027 (S05E48 L=102)   603    145        2 ( 257/1002- 364/1054)   6.3 E52-E43 S05-S02
      9 1577  776 (S17W37 L=187)   957    291        1 (1461/ 656-1654/ 895)  53.7 W30-W41 S24-S08
     10 1399  852 (S13W24 L=173)   511     78        1 (1342/ 821-1439/ 888)  33.7 W20-W26 S14-S10
     11  288  706 (S24E53 L= 97)   989    225        0 ( 241/ 647- 344/ 751)  26.6 E58-E46 S27-S21
     12  590  568 (S33E30 L=119)   745    221        2 ( 494/ 537- 669/ 605)  11.3 E38-E23 S35-S30
  

    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