Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  227 1383 (N25E58 L= 68)   499    110        1 ( 196/1367- 258/1395)  -7.1 E62-E54 N23-N26
      2  513 1284 (N16E32 L= 94)   350     91        0 ( 466/1260- 545/1302) -18.4 E35-E29 N14-N17
      3  168 1284 (N19E64 L= 63)   491    123        0 ( 141/1268- 196/1298) -15.9 E67-E60 N17-N20
      4  237 1224 (N15E55 L= 71)   452    111        0 ( 202/1207- 278/1239) -14.0 E58-E50 N13-N15
      5 1200 1158 (N05W11 L=137)   276     77        0 (1162/1117-1230/1191) -53.1 W09-W12 N02-N06
      6 1506 1017 (S05W30 L=156)   244     62        0 (1487/ 981-1525/1050) -63.4 W28-W31 S06-S02
      7 1682 1001 (S06W43 L=169)   422     94        1 (1651/ 955-1706/1042) -59.0 W40-W45 S08-S03
      8  126  940 (S02E70 L= 56)   567    173        0 (  92/ 922- 158/ 963) -10.3 E76-E64 S02-N00
      9 1846  678 (S26W62 L=188)   806     77        0 (1820/ 643-1875/ 717)  59.0 W59-W64 S28-S23
     10  404  689 (S19E44 L= 82)   293     94        0 ( 367/ 657- 426/ 712) -33.7 E47-E42 S20-S17
     11  376  607 (S24E49 L= 77)   433    100        2 ( 352/ 577- 416/ 653) -68.2 E51-E46 S25-S20
     12 1078  489 (S36W01 L=127)   621    142        1 (1011/ 467-1155/ 508)   5.7 E04-W06 S37-S34
  

    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