Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1137 1626 (N36W24 L=177)  1299    225        3 (1043/1570-1252/1658) -14.9 W16-W31 N32-N37
      2 1776 1380 (N05W60 L=213)   861     93        1 (1735/1338-1815/1424) -63.4 W56-W63 N02-N08
      3  161 1306 (N37E60 L= 93)   818    280        1 ( 121/1268- 202/1323) -20.0 E65-E54 N34-N38
      4 1764 1179 (S05W51 L=204)   461    107        0 (1735/1149-1803/1205)   8.1 W48-W55 S06-S03
      5  450 1142 (N21E33 L=120)   281     63        1 ( 422/1121- 482/1158) -26.6 E35-E31 N20-N22
      6 1318  979 (S06W16 L=169)  1099    344        3 (1237/ 800-1378/1143)  45.0 W07-W22 S14-N01
      7 1786 1091 (S11W51 L=204)   345     78        0 (1750/1066-1812/1110) -31.0 W48-W53 S12-S09
      8 1777  950 (S19W47 L=201)   412    117        3 (1730/ 913-1809/ 994)  23.2 W42-W50 S20-S16
      9  491  955 (N10E33 L=120)   232     66        1 ( 474/ 931- 508/ 987) -90.0 E34-E32 N08-N11
     10  925  447 (S29E21 L=132)   875    213        1 ( 886/ 357- 949/ 529)  54.0 E26-E17 S34-S23
     11  534  399 (S24E49 L=104)  1576    404        8 ( 484/ 286- 570/ 512)  58.0 E56-E41 S32-S16
     12 1346  474 (S38W08 L=161)   522    100        0 (1297/ 426-1388/ 512)  26.6 W03-W11 S40-S35
  

    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