Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  507 1720 (N70E37 L=115)   790    100        3 ( 474/1680- 549/1751)  90.0 E45-E26 N66-N73
      2 1387 1612 (N30W44 L=196)   357     39        1 (1358/1591-1412/1634) -11.3 W41-W46 N28-N30
      3 1586 1600 (N22W59 L=210)  1073    176        1 (1500/1579-1671/1623) -26.6 W52-W67 N18-N25
      4 1622 1412 (N11W50 L=201)   279     70        0 (1592/1387-1651/1437) -63.4 W48-W51 N09-N13
      5  523 1274 (N34E27 L=124)   283    101        0 ( 475/1253- 566/1296) -32.0 E31-E23 N31-N36
      6  492 1098 (N24E32 L=119)   457    148        0 ( 436/1065- 558/1146) -41.2 E35-E27 N20-N27
      7  252  838 (N13E58 L= 93)   615    106        1 ( 210/ 819- 307/ 867) -16.7 E63-E52 N10-N14
      8  425  769 (N06E44 L=107)   349     88        1 ( 402/ 722- 444/ 805) -90.0 E45-E42 N02-N08
      9  559  548 (S10E42 L=110)   294     82        1 ( 528/ 526- 600/ 569)   0.0 E45-E38 S10-S08
     10  413  471 (S12E58 L= 93)   648     81        0 ( 394/ 429- 436/ 510)  32.0 E61-E53 S15-S09
     11 1190  427 (S34E07 L=145)   299     60        2 (1153/ 398-1216/ 449) -26.6 E08-E04 S35-S31
  

    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