Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  969 1709 (N47E15 L=289)   366     73        1 ( 936/1683-1000/1734)  59.0 E17-E12 N44-N49
      2  321 1476 (N22E63 L=242)  2197    202        5 ( 282/1408- 357/1580) -59.7 E69-E56 N16-N28
      3 1341 1324 (N23W18 L=322)   211     92        0 (1297/1306-1374/1343)  -9.5 W14-W20 N22-N24
      4  489 1307 (N13E39 L=265)   331     92        1 ( 447/1293- 536/1320)  -8.1 E42-E35 N12-N13
      5  617  897 (S11E25 L=280)   470    114        2 ( 576/ 855- 673/ 928)  39.8 E27-E20 S13-S08
      6  400  841 (S17E41 L=263)   330     40        2 ( 376/ 815- 420/ 873) -90.0 E43-E39 S18-S14
      7  187  824 (S20E66 L=238)   724     68        1 ( 158/ 792- 203/ 854)  24.0 E72-E62 S22-S18
      8  645  794 (S17E22 L=282)   420    151        2 ( 579/ 772- 703/ 833) -14.0 E26-E18 S18-S14
      9 1341  726 (S14W23 L=328)   539    176        0 (1321/ 659-1371/ 798) -76.0 W21-W26 S17-S09
     10  589  641 (S28E27 L=278)   276    112        0 ( 574/ 594- 602/ 699)  74.1 E27-E25 S30-S23
     11 1341  386 (S37W34 L=338)   781     97        3 (1295/ 353-1381/ 417)  15.9 W29-W37 S39-S34
     12 1269  235 (S52W40 L=344)   888    103        0 (1242/ 187-1297/ 285) -50.7 W35-W46 S57-S46
  

    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