Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  357 1465 (N50E43 L=  8)   635    110        1 ( 302/1441- 414/1481) -13.2 E51-E34 N47-N51
      2  734 1337 (N31E09 L= 42)  1069    234        1 ( 620/1272- 834/1369) -31.6 E16-E03 N25-N35
      3  523 1305 (N34E25 L= 26)   355     70        0 ( 486/1279- 549/1331)  18.4 E28-E22 N33-N35
      4  497 1148 (N25E29 L= 21)   294     62        0 ( 463/1126- 531/1168) -36.9 E31-E27 N23-N27
      5  870 1099 (N13E07 L= 44)   214    117        0 ( 848/1055- 909/1148)  38.7 E08-E03 N10-N15
      6  685 1066 (N16E18 L= 32)   341    125        1 ( 621/1046- 747/1088) -35.5 E22-E14 N13-N18
      7  280 1031 (N23E50 L=  1)   475     70        0 ( 253/ 999- 314/1058) -45.0 E52-E47 N20-N25
      8  749  617 (S10E27 L= 24)   225     88        0 ( 706/ 596- 795/ 633)  -8.1 E30-E23 S11-S09
      9  616  340 (S24E50 L=  1)  1792    333        4 ( 514/ 230- 669/ 495)  64.5 E58-E45 S33-S12
     10 1035  316 (S38E23 L= 28)  2377    423        4 ( 867/ 223-1171/ 437)   8.9 E40-E08 S41-S32
  

    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