Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1625 1481 (N13W53 L=128)  1394    207        5 (1560/1391-1680/1551)  85.2 W49-W55 N06-N18
      2  253 1462 (N52E63 L= 13)  1304    105        2 ( 202/1436- 300/1484)   0.0 E79-E53 N49-N53
      3  835 1427 (N33E00 L= 76)   328     67        1 ( 802/1399- 867/1445) -45.0 E01-W02 N31-N34
      4  967 1386 (N27W07 L= 82)   641    140        2 ( 916/1352-1040/1408) -14.0 W03-W11 N25-N29
      5 1283 1354 (N16W24 L=100)   648    339        0 (1147/1334-1448/1378) -24.0 W16-W34 N12-N20
      6  286  944 (N18E51 L= 25)  1604    370        3 ( 260/ 800- 334/1121)  81.5 E54-E46 N08-N28
      7 1421 1019 (S06W22 L= 98)   273     57        1 (1400/ 990-1438/1056)  63.4 W20-W23 S07-S03
      8 1908  999 (S22W64 L=140)   566     39        1 (1886/ 972-1930/1020)  18.4 W60-W68 S23-S20
      9 1715  850 (S24W40 L=115)   957    227        1 (1675/ 733-1763/ 940)  54.0 W34-W45 S30-S19
     10  344  459 (S12E67 L=  8)   625     70        2 ( 322/ 422- 372/ 484)  26.6 E72-E62 S14-S10
     11  550  329 (S24E57 L= 19)  1117    176        0 ( 511/ 272- 571/ 419)  25.3 E68-E49 S27-S18
     12  807  282 (S34E40 L= 35)   349     32        0 ( 791/ 262- 827/ 309)  53.1 E42-E38 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