Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1448 1610 (N45W29 L=119)   445     71        0 (1412/1594-1486/1627)  -7.1 W25-W33 N43-N46
      2  387 1557 (N27E60 L= 30)   977    115        3 ( 364/1491- 413/1613) -48.8 E65-E56 N23-N31
      3 1221 1410 (N27W09 L= 99)   926    290        0 (1116/1304-1277/1488) -47.7 W01-W13 N20-N31
      4  322 1449 (N20E60 L= 30)   566     77        1 ( 290/1425- 352/1466)  -6.3 E65-E56 N18-N20
      5 1699 1392 (N31W50 L=140)   369     42        1 (1679/1364-1718/1419) -56.3 W47-W51 N29-N33
      6 1347 1236 (N17W19 L=109)   352    101        1 (1289/1213-1388/1263)  -8.1 W14-W21 N16-N18
      7  293  735 (S25E53 L= 37)   431     78        1 ( 254/ 715- 339/ 747)   5.7 E57-E47 S26-S23
      8 1084  702 (S18W07 L= 97)   432     88        2 (1043/ 669-1129/ 739)  45.0 W04-W09 S20-S15
      9 1529  601 (S20W40 L=131)   643    141        1 (1476/ 570-1584/ 627)  20.6 W36-W45 S22-S18
     10  359  544 (S38E52 L= 38)  1219    159        1 ( 299/ 502- 429/ 586)  16.9 E64-E41 S42-S33
     11 1077  187 (S60W24 L=114)   688     97        0 (1055/ 142-1092/ 246) -47.3 W18-W32 S66-S53
  

    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