Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1005 1620 (N31E18 L=232)  2015    352        7 ( 892/1516-1092/1688) -17.5 E28-E09 N25-N34
      2  791 1576 (N22E30 L=221)  1059    197        3 ( 730/1498- 864/1655)  85.6 E31-E27 N16-N29
      3 1500 1325 (N26W22 L=272)  2403    536        5 (1365/1094-1566/1495) -49.6 W09-W29 N13-N35
      4  188 1226 (S11E65 L=185)   853     79        1 ( 160/1194- 216/1269)  56.3 E68-E61 S13-S07
      5  152 1105 (S18E68 L=183)  2062    239        2 (  83/1026- 233/1193)  25.6 E84-E58 S23-S11
      6  410  978 (S21E38 L=213)   697    189        1 ( 344/ 920- 503/1018)  39.3 E42-E31 S25-S16
      7 1776  957 (N12W50 L=301)   407     84        1 (1759/ 913-1788/1007) -80.5 W48-W51 N09-N15
      8  852  852 (S18E06 L=245)   355    151        0 ( 821/ 796- 902/ 895)  71.6 E07-E03 S20-S14
      9  425  856 (S28E35 L=216)   444     93        1 ( 378/ 843- 482/ 874)  12.5 E39-E30 S28-S26
     10 1310  788 (S10W22 L=273)   365    127        0 (1251/ 753-1366/ 810)  14.0 W18-W26 S11-S08
     11  744  392 (S51E00 L=250)  1688    274        3 ( 640/ 311- 839/ 490) -17.7 E11-W11 S54-S46
     12  946  245 (S55W27 L=277)   386     34        1 ( 923/ 225- 972/ 262) -78.7 W23-W29 S57-S52
  

    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