Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1598 1641 (N36W62 L=169)   910     94        0 (1574/1600-1635/1670)   7.1 W55-W71 N34-N37
      2 1740 1377 (N17W59 L=165)  1238    275        2 (1689/1275-1793/1470) -85.6 W55-W64 N10-N23
      3 1365 1444 (N28W29 L=136)   237     75        0 (1330/1420-1396/1469) -53.1 W27-W31 N25-N29
      4 1704 1329 (N15W53 L=160)   806    140        3 (1662/1297-1755/1386) -63.4 W48-W57 N12-N19
      5  312 1278 (N27E53 L= 54)  1224    279        4 ( 222/1240- 433/1311)  -4.8 E64-E40 N24-N28
      6  720 1192 (N19E18 L= 88)   253    139        0 ( 705/1135- 754/1253)  66.8 E19-E15 N15-N22
      7  234 1176 (N20E60 L= 46)   821    126        1 ( 209/1123- 261/1233)  60.3 E62-E57 N16-N23
      8  773  816 (S05E18 L= 89)   697    224        1 ( 749/ 720- 807/ 923)  84.8 E19-E16 S10-N01
      9 1749  856 (S14W50 L=156)   508     23        0 (1723/ 833-1771/ 879)  56.3 W47-W51 S15-S12
     10  429  743 (S07E43 L= 63)   326    133        0 ( 393/ 695- 465/ 806) -69.4 E45-E41 S10-S02
     11  409  544 (S20E52 L= 55)   356     95        0 ( 394/ 496- 421/ 573)  59.0 E54-E49 S23-S18
     12 1360  487 (S34W18 L=125)   357     73        0 (1315/ 468-1398/ 504) -23.2 W14-W21 S35-S31
  

    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