Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  681 1724 (N41E43 L= 89)   425     72        1 ( 643/1704- 717/1741)  -6.3 E46-E37 N39-N42
      2  843 1657 (N38E25 L=106)   417     98        1 ( 808/1633- 898/1677)  26.6 E28-E20 N36-N39
      3 1304 1428 (N29W13 L=144)  2193    463        1 (1098/1327-1490/1502)  -8.1 E01-W27 N24-N31
      4 1610 1226 (N20W37 L=168)   322    107        0 (1568/1184-1640/1255) -33.7 W33-W39 N17-N21
      5  491 1126 (S01E35 L= 96)   869    178        1 ( 447/1067- 531/1206)  69.4 E37-E32 S04-N04
      6 1500 1172 (N15W28 L=160)   359    113        1 (1475/1123-1542/1200)  45.0 W26-W31 N12-N17
      7 1795  966 (N07W55 L=186)   424    117        0 (1774/ 903-1813/1016) -56.3 W52-W57 N02-N09
      8  775  785 (S17E12 L=119)   385    136        0 ( 722/ 754- 809/ 839)  45.0 E15-E10 S19-S13
      9 1560  712 (S11W40 L=171)   776    202        0 (1469/ 658-1627/ 777)  38.7 W33-W43 S15-S06
     10 1594  541 (S22W49 L=180)  2188    475        0 (1396/ 488-1732/ 670)  36.0 W34-W57 S27-S11
  

    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