Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1157 1730 (N52W20 L= 51)   330     60        1 (1126/1712-1195/1741)   0.0 W16-W23 N50-N53
      2 1579 1546 (N34W50 L= 82)  1351    246        2 (1442/1520-1677/1589) -17.7 W38-W61 N30-N38
      3  520 1427 (N32E37 L=355)   836    182        1 ( 472/1371- 572/1501) -36.9 E43-E31 N28-N37
      4 1362 1264 (N17W24 L= 55)   367    153        0 (1302/1224-1429/1308) -40.6 W19-W27 N14-N20
      5  512 1199 (N17E34 L=358)   393    143        0 ( 478/1143- 559/1245)  45.0 E36-E30 N13-N20
      6  697 1105 (N11E20 L= 11)   369    187        1 ( 634/1032- 744/1164) -48.8 E24-E17 N06-N14
      7 1012 1104 (N10E00 L= 31)   481     96        1 ( 975/1058-1043/1152) -80.5 E02-W01 N06-N12
      8 1659 1095 (N04W43 L= 75)   304     63        0 (1629/1076-1690/1114)   0.0 W40-W46 N03-N05
      9 1364  783 (S12W20 L= 52)   439    124        2 (1324/ 741-1422/ 820)  33.7 W17-W24 S14-S10
     10  372  722 (S12E48 L=344)  1120    191        1 ( 342/ 633- 435/ 780)  48.0 E52-E41 S18-S08
  

    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