Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  319 1462 (N21E62 L=241)  1644    140        3 ( 293/1406- 358/1530) -66.0 E68-E55 N16-N25
      2 1349 1323 (N23W19 L=322)   244     87        0 (1300/1304-1386/1339)  -9.5 W15-W21 N22-N23
      3  501 1306 (N13E38 L=265)   315     89        0 ( 460/1290- 547/1317)  -8.1 E41-E34 N12-N13
      4  647  901 (S10E23 L=281)   353    106        2 ( 613/ 872- 690/ 919)  21.8 E24-E19 S12-S08
      5  665  789 (S17E21 L=282)   420    158        1 ( 598/ 769- 719/ 829) -14.0 E25-E17 S18-S14
      6 1354  723 (S14W24 L=328)   528    160        1 (1332/ 660-1387/ 792) -82.9 W22-W27 S17-S09
      7 1084  718 (S17W07 L=310)   368    164        0 (1031/ 691-1139/ 758)  33.7 W03-W09 S18-S13
      8  602  644 (S27E26 L=278)   262     84        1 ( 584/ 604- 616/ 687)  71.6 E27-E24 S30-S24
      9 1711  641 (S16W56 L=359)   396     57        1 (1684/ 619-1732/ 660)  14.0 W52-W58 S17-S14
     10 1351  383 (S37W35 L=338)   746     88        2 (1309/ 345-1389/ 412)  18.4 W30-W38 S40-S34
     11 1271  234 (S52W40 L=343)   823    103        1 (1236/ 187-1303/ 281) -50.7 W36-W47 S57-S46
  

    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