Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1182 1707 (N47E06 L=254)   452    100        0 (1128/1695-1232/1722)  26.6 E10-E01 N44-N48
      2  746 1590 (N27E32 L=228)  2275    357        9 ( 656/1438- 812/1719) -34.8 E46-E23 N18-N34
      3  564 1382 (N11E37 L=224)   299     99        2 ( 536/1347- 608/1403)  53.1 E38-E34 N08-N12
      4 1688 1280 (N28W42 L=302)   755    170        2 (1620/1222-1769/1321)  -3.8 W35-W50 N25-N29
      5 1250 1212 (N14W10 L=270)   552    161        4 (1219/1138-1276/1261) -45.0 W07-W12 N10-N17
      6 1027  861 (S11W03 L=264)   495    123        1 ( 959/ 839-1074/ 888)  26.6 W00-W06 S12-S08
      7 1074  697 (S20W10 L=270)   526    219        1 ( 993/ 650-1157/ 728)  33.7 W05-W14 S22-S16
      8  443  637 (S38E35 L=225)   287     40        1 ( 422/ 615- 457/ 665) -90.0 E37-E33 S39-S35
      9 1258  497 (S28W28 L=289)   330    102        0 (1227/ 462-1283/ 540)  90.0 W27-W29 S31-S25
     10  468  390 (S58E41 L=220)  1113    181        1 ( 398/ 368- 546/ 414)  10.1 E57-E27 S61-S53
     11 1008  376 (S43W16 L=277)   482    141        0 ( 945/ 356-1066/ 411)  49.4 W11-W18 S45-S38
     12 1328  355 (S36W42 L=302)   813    144        1 (1264/ 333-1390/ 385)  32.0 W36-W45 S38-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