Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  303 1559 (N09E68 L=258)  1776    320        4 ( 259/1522- 358/1616)   8.7 E75-E61 N06-N12
      2  666 1545 (N14E36 L=291)   304     79        2 ( 642/1514- 694/1574) -26.6 E38-E32 N12-N15
      3  749 1394 (N07E25 L=301)   204     45        1 ( 734/1371- 775/1417)  90.0 E26-E23 N05-N08
      4  139 1160 (S17E68 L=259)  1343    284        3 (  86/1115- 189/1194)   8.7 E74-E61 S20-S14
      5  640 1151 (S09E25 L=301)   276     63        1 ( 619/1126- 660/1189) -36.9 E27-E23 S10-S07
      6  367 1072 (S19E43 L=284)   941    273        1 ( 266/1011- 456/1141)  -9.5 E52-E34 S21-S16
      7  195 1015 (S25E57 L=269)   840    217        1 ( 148/ 948- 243/1082)  52.1 E61-E53 S29-S20
      8  532 1033 (S18E30 L=297)   263     73        0 ( 495/1010- 565/1063)  51.3 E31-E27 S20-S15
      9  595  967 (S21E24 L=302)   657    148        0 ( 518/ 942- 657/ 996)  16.7 E29-E19 S22-S18
     10  458  766 (S36E31 L=295)   290    113        1 ( 412/ 745- 505/ 791)   0.0 E35-E27 S36-S34
     11 1356  736 (S14W27 L=353)   217     84        1 (1340/ 695-1379/ 775) -53.1 W25-W28 S15-S11
     12 1396  399 (S31W45 L= 12)  2586    583        4 (1167/ 349-1603/ 456)  17.9 W27-W61 S35-S24
  

    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