Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1550 1363 (N11W40 L=313)   671    137        1 (1517/1295-1598/1433) -90.0 W38-W41 N06-N15
      2  390 1255 (N27E38 L=235)   277     34        0 ( 378/1226- 403/1282)  71.6 E39-E36 N25-N28
      3  308 1228 (N27E45 L=228)   336     79        1 ( 272/1204- 343/1250)   9.5 E48-E41 N25-N27
      4  766  953 (N02E16 L=257)   265     59        0 ( 747/ 925- 790/ 989)  76.0 E17-E15 N00-N04
      5 1890  823 (S28W60 L=333)  1393    189        2 (1842/ 768-1935/ 869)  21.8 W54-W65 S30-S25
      6  650  810 (S04E26 L=247)   416    127        0 ( 587/ 785- 692/ 852)   7.1 E30-E22 S04-S01
      7  922  357 (S37E23 L=250)   321     44        1 ( 898/ 340- 952/ 379)   0.0 E25-E20 S38-S35
      8  689  253 (S40E47 L=227)   373     23        0 ( 671/ 234- 706/ 271)  11.3 E48-E43 S41-S38
  

    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