Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1601 1517 (N37W46 L=203)  1407    178        3 (1522/1477-1676/1546)  -7.1 W37-W54 N34-N39
      2 1347 1450 (N31W21 L=178)   326     83        0 (1313/1423-1383/1481) -45.0 W19-W23 N29-N33
      3  379 1444 (N25E53 L=104)   403     58        1 ( 353/1422- 410/1465)  31.0 E55-E50 N23-N26
      4 1424 1373 (N26W27 L=183)   359     62        1 (1391/1354-1457/1395)   0.0 W23-W29 N25-N27
      5 1809 1329 (N24W64 L=221)  1491    231        1 (1724/1262-1878/1391) -27.9 W54-W72 N19-N28
      6 1748  855 (S05W52 L=209)   511     21        0 (1723/ 833-1775/ 877)  26.6 W50-W55 S06-S04
      7  292  737 (S20E54 L=103)   677     99        1 ( 249/ 709- 335/ 766)  -6.3 E58-E49 S21-S18
      8 1665  683 (S17W48 L=205)   378     77        0 (1647/ 644-1683/ 722)  90.0 W46-W49 S19-S13
      9  274  675 (S24E59 L= 98)   440     62        0 ( 250/ 653- 305/ 696)  18.4 E63-E54 S26-S22
     10 1780  605 (S22W66 L=223)  1288    171        2 (1747/ 542-1808/ 686) -83.7 W61-W73 S26-S16
     11  375  585 (S30E50 L=107)  1331    217        0 ( 308/ 495- 418/ 684) -69.0 E54-E46 S36-S23
     12  624  507 (S34E28 L=129)   441     58        1 ( 592/ 492- 668/ 527) -11.3 E30-E24 S34-S32
  

    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