Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  421 1514 (N06E53 L=206)   333     58        0 ( 403/1483- 446/1541)  90.0 E54-E51 N04-N08
      2 1542 1507 (N38W19 L=278)   562     63        0 (1505/1483-1587/1531)   0.0 W15-W23 N36-N39
      3  236 1285 (S10E61 L=198)  1012    145        1 ( 168/1248- 275/1336)  24.4 E69-E57 S12-S06
      4  488 1240 (S08E38 L=221)   619    168        0 ( 460/1163- 537/1308)  84.3 E39-E35 S12-S02
      5 1173 1249 (N10W02 L=261)   249     63        0 (1145/1219-1195/1281) -45.0 W00-W03 N08-N11
      6  330 1022 (S24E45 L=214)   855    201        0 ( 289/ 936- 406/1082)  54.0 E48-E39 S29-S18
      7 1789  818 (N06W55 L=314)   636    114        1 (1760/ 759-1815/ 871) -81.9 W52-W56 N02-N09
      8 1630  767 (S02W43 L=302)   425    113        0 (1580/ 731-1673/ 807)  56.3 W40-W45 S05-N01
      9 1250  753 (S14W20 L=279)   270     49        1 (1224/ 729-1278/ 777)  56.3 W18-W21 S15-S12
     10  434  584 (S49E32 L=227)  1051    218        2 ( 337/ 557- 526/ 635)   2.6 E44-E21 S51-S45
     11 1350  479 (S26W38 L=297)   550    124        0 (1284/ 458-1406/ 499)  11.3 W33-W43 S27-S24
     12  994  329 (S47W26 L=285)  1325    260        0 ( 879/ 268-1103/ 382)   4.8 W14-W38 S48-S45
  

    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