Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  763 1651 (N52E09 L=249)   411     68        2 ( 733/1613- 794/1675)  45.0 E12-E06 N49-N54
      2 1602 1549 (N26W56 L=314)   349     63        1 (1576/1519-1630/1572) -63.4 W54-W57 N24-N28
      3 1748 1334 (N11W58 L=316)  1090    225        1 (1687/1271-1815/1419) -61.4 W54-W63 N05-N16
      4  194 1076 (N19E64 L=195)  1046    175        3 ( 146/1048- 271/1111)   6.3 E71-E53 N17-N21
      5 1173  894 (S04W07 L=265)   384    110        2 (1117/ 875-1227/ 906) -15.9 W03-W10 S05-S02
      6  273  841 (N04E56 L=202)   325     96        0 ( 255/ 802- 293/ 882)  38.7 E58-E53 N02-N06
      7 1299  504 (S31W10 L=268)  3063    647        1 (1025/ 363-1488/ 692)  19.9 E12-W24 S37-S21
      8  274  616 (S10E64 L=194)   462     45        1 ( 257/ 594- 297/ 640) -45.0 E66-E61 S11-S08
      9  529  546 (S16E42 L=217)   326    101        0 ( 506/ 504- 563/ 580)  29.7 E45-E38 S19-S14
     10  747  400 (S29E32 L=227)   550    131        1 ( 731/ 341- 776/ 466)  68.2 E33-E28 S34-S24
     11  538  359 (S30E52 L=206)   810    173        0 ( 514/ 294- 584/ 437)  30.1 E61-E42 S35-S24
     12  813  266 (S41E34 L=224)   401     63        1 ( 783/ 242- 845/ 290)   0.0 E37-E30 S43-S39
  

    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