Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1593 1531 (N33W48 L= 40)   540    102        1 (1546/1502-1641/1560) -33.7 W43-W51 N30-N35
      2  746 1341 (N25E18 L=334)  3682   1129       10 ( 336/1101-1062/1552) -27.0 E57-W02 N08-N39
      3 1055 1465 (N32W04 L=356)   225     53        2 (1023/1448-1087/1481)  11.3 W01-W06 N31-N33
      4 1453 1420 (N27W32 L= 24)   539    173        0 (1375/1388-1534/1436) -18.4 W26-W38 N24-N28
      5 1688 1233 (N13W47 L= 40)   913    237        0 (1609/1150-1773/1302) -42.5 W41-W54 N07-N18
      6  341 1176 (N15E48 L=304)   764    200        2 ( 323/1092- 367/1280) -85.2 E50-E45 N09-N21
      7  164  925 (S01E68 L=284)   745    197        0 ( 128/ 905- 210/ 943)   8.1 E75-E61 S03-N00
      8 1420  858 (S08W25 L= 17)   509    127        2 (1367/ 820-1477/ 906) -40.6 W21-W28 S10-S04
      9  997  771 (S11E03 L=350)   352     78        2 ( 974/ 731-1015/ 815)  78.7 E04-E01 S13-S08
     10 1562  761 (S15W35 L= 27)   315     49        1 (1531/ 743-1585/ 783)  26.6 W32-W36 S15-S13
     11  879  642 (S19E11 L=341)   630    210        0 ( 794/ 575- 954/ 720)  36.9 E17-E05 S23-S14
     12 1097  502 (S30W03 L=355)  1818    388        4 ( 967/ 385-1274/ 604)  32.5 E07-W15 S38-S23
  

    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