Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1415 1706 (N31W52 L= 80)   681    106        1 (1366/1681-1468/1741) -63.4 W50-W55 N28-N34
      2  542 1559 (N51E21 L=  7)   365     39        0 ( 517/1542- 571/1576) -18.4 E23-E17 N49-N53
      3  202 1313 (N43E62 L=325)   508     49        0 ( 168/1299- 234/1326)   0.0 E69-E56 N41-N43
      4 1217 1231 (N11W16 L= 44)   187    109        0 (1190/1188-1243/1279) -90.0 W15-W17 N08-N14
      5  531 1162 (N25E27 L=  1)   295     89        0 ( 492/1138- 575/1177) -33.7 E29-E23 N22-N26
      6  176 1113 (N30E62 L=325)  1089    150        1 ( 121/1051- 229/1162) -26.6 E73-E55 N25-N33
      7  337  997 (N20E45 L=342)   377     82        1 ( 295/ 980- 379/1012)  -7.1 E49-E41 N18-N20
      8 1425  510 (S37W11 L= 38)   882    154        1 (1370/ 451-1470/ 583)  41.2 W05-W15 S40-S32
      9  464  482 (S12E53 L=335)   443     79        2 ( 436/ 449- 493/ 510) -76.0 E54-E50 S14-S10
     10  874  343 (S31E31 L=356)  1888    366        6 ( 737/ 219- 971/ 465)  19.0 E51-E19 S36-S24
     11 1187  368 (S39E10 L= 17)  1208    284        2 (1067/ 309-1309/ 431)   2.9 E20-W00 S41-S37
  

    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