Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  894 1713 (N31E32 L= 31)   405     57        1 ( 862/1686- 919/1743) -80.5 E33-E29 N27-N33
      2  378 1461 (N02E55 L=  7)   497     83        2 ( 352/1427- 421/1487)  53.1 E57-E52 N00-N04
      3  759 1329 (N03E23 L= 39)   428    141        0 ( 728/1278- 800/1397) -32.0 E27-E19 N01-N06
      4 1011 1073 (S04E02 L= 60)   593    204        0 ( 942/ 986-1065/1152)  84.8 E03-W00 S10-N01
      5  273  956 (S29E51 L= 12)   727     99        1 ( 232/ 918- 317/ 989)  29.1 E55-E46 S32-S26
      6 1022  815 (S18W06 L= 68)   702    253        0 ( 940/ 770-1126/ 899)  59.7 W02-W09 S22-S10
      7  750  832 (S25E11 L= 52)   501    232        1 ( 655/ 789- 856/ 854)  12.1 E16-E02 S27-S22
      8  506  768 (S35E27 L= 36)   283     85        0 ( 461/ 751- 550/ 782)  29.7 E30-E23 S37-S33
      9  606  635 (S41E16 L= 47)  1236    320        0 ( 482/ 569- 732/ 720)  41.6 E26-E08 S47-S31
     10 1576  366 (S24W66 L=129)  1741    134        3 (1515/ 340-1644/ 396)  12.5 W57-W75 S26-S20
     11 1391  293 (S34W58 L=120)   517     71        2 (1364/ 257-1410/ 329) -16.7 W53-W63 S36-S32
     12 1234  205 (S46W58 L=121)  2040    160        2 (1193/ 152-1298/ 276)  -6.3 W47-W74 S48-S41
  

    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