Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1332 1603 (N36W04 L= 23)   455     79        2 (1305/1572-1362/1635) -21.8 W00-W05 N33-N38
      2  240 1482 (N04E71 L=308)   909    141        3 ( 213/1453- 279/1503)  38.7 E76-E66 N01-N05
      3  162 1189 (S15E66 L=314)  1476    394        3 (  86/1099- 238/1290)  47.5 E73-E60 S20-S08
      4  341 1108 (S17E46 L=334)  1347    294        1 ( 276/ 989- 386/1254)  76.8 E49-E43 S24-S07
      5  617 1107 (S12E26 L=354)   561    107        2 ( 580/1046- 660/1137) -38.7 E28-E23 S14-S09
      6  763 1079 (S10E16 L=  4)   188     73        0 ( 723/1060- 799/1095)   0.0 E18-E13 S10-S09
      7 1907  996 (N17W62 L= 82)  1454    209        5 (1881/ 946-1943/1043) -81.9 W58-W67 N13-N20
      8 1694  767 (S03W49 L= 68)  1356    261        1 (1610/ 719-1778/ 852)  -3.2 W39-W57 S04-N00
      9 1110  761 (S20W12 L= 32)   483    158        0 (1085/ 686-1133/ 841) -66.0 W08-W14 S24-S15
     10  772  711 (S31E07 L= 12)  1006    339        1 ( 596/ 670- 919/ 741)  12.8 E19-W03 S35-S28
     11 1509  425 (S26W51 L= 71)  2554    475        5 (1358/ 351-1716/ 517)  45.0 W42-W65 S34-S15
     12 1630  386 (S24W64 L= 84)  1562    214        1 (1582/ 335-1663/ 453) -58.0 W60-W66 S27-S19
  

    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