Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  541 1490 (N09E41 L=280)  2758    553        3 ( 487/1259- 619/1719) -63.4 E47-E32 S03-N25
      2  240 1485 (N05E71 L=251)   659    160        1 ( 213/1468- 263/1503) -10.3 E77-E66 N04-N06
      3 1787 1461 (N42W46 L=  8)   488     86        0 (1768/1433-1811/1482)  45.0 W43-W49 N39-N44
      4 1665 1384 (N31W33 L=355)   379     68        1 (1625/1356-1695/1404) -11.3 W29-W35 N29-N32
      5 1338 1294 (N16W12 L=333)   479    146        0 (1304/1220-1369/1371) -53.1 W08-W14 N11-N19
      6  181 1232 (S11E65 L=256)  1467    275        1 ( 139/1179- 217/1292)  54.5 E71-E61 S14-S06
      7  248 1166 (S14E56 L=266)  1197    298        0 ( 162/1142- 333/1185)   6.3 E65-E47 S15-S11
      8 1378 1084 (N05W19 L=340)   287     84        1 (1344/1061-1416/1103)  14.0 W16-W20 N02-N05
      9 1742  967 (N08W46 L=  8)   500     81        1 (1724/ 923-1768/1012) -45.0 W43-W49 N05-N10
     10  751  918 (S19E14 L=308)   308    109        0 ( 705/ 900- 820/ 934)  29.7 E16-E09 S21-S16
     11 1451  688 (S15W34 L=356)   600    132        0 (1379/ 640-1517/ 740)  60.9 W31-W37 S19-S10
     12 1699  557 (S14W59 L= 21)   732    132        0 (1649/ 506-1735/ 612)  76.0 W56-W61 S18-S10
     13 1582  428 (S25W56 L= 18)  1626    215        3 (1527/ 391-1645/ 477)  32.0 W50-W61 S27-S21
  

    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