Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1034 1817 (N43E31 L=199)   456     96        0 (1000/1769-1058/1859) -26.6 E37-E25 N39-N46
      2  971 1457 (N16E16 L=214)   799    225        0 ( 904/1358-1044/1521) -26.6 E21-E08 N12-N19
      3 1840 1062 (N22W54 L=283)  1070    186        0 (1778/ 980-1872/1118) -24.4 W46-W59 N17-N24
      4  378 1025 (S23E42 L=188)   575    174        0 ( 326/ 965- 442/1072)  45.0 E45-E37 S27-S18
      5  593 1026 (S18E25 L=204)   617    243        0 ( 491/ 979- 710/1066)  37.6 E31-E18 S22-S12
      6  718  864 (S24E13 L=216)   469    176        1 ( 648/ 828- 803/ 903)  45.0 E17-E09 S27-S19
      7  257  849 (S36E53 L=177)   361     69        0 ( 241/ 813- 277/ 886)  59.0 E54-E49 S38-S33
      8  897  675 (S29W03 L=232)  2097    691        4 ( 653/ 556-1222/ 748)   9.7 E14-W27 S34-S24
      9 1439  565 (S19W40 L=269)   457    115        1 (1390/ 530-1491/ 601)  63.4 W37-W42 S21-S15
     10  323  488 (S61E58 L=171)   598     38        1 ( 293/ 472- 347/ 504)   3.8 E68-E50 S62-S59
     11 1496  327 (S29W63 L=292)   981    141        2 (1436/ 290-1555/ 376)  82.9 W59-W65 S32-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