Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  453 1269 (N35E33 L=159)   643    174        2 ( 383/1229- 527/1308) -28.3 E40-E27 N30-N38
      2  822 1290 (N27E05 L=187)   387    125        1 ( 755/1267- 883/1306) -12.5 E10-E01 N25-N28
      3  283 1228 (N35E51 L=141)  1307    248        2 ( 209/1155- 371/1298)  23.6 E60-E41 N31-N38
      4 1250 1252 (N14W19 L=211)   335     93        1 (1203/1220-1288/1282) -51.3 W16-W20 N11-N16
      5 1652 1190 (S01W43 L=235)   364     84        1 (1602/1163-1682/1207) -18.4 W39-W45 S02-N01
      6 1236 1077 (N05W13 L=206)   286    119        0 (1206/1033-1272/1122)  26.6 W10-W16 N02-N06
      7 1782  923 (S20W48 L=240)   895    192        1 (1715/ 842-1839/ 988)  15.9 W40-W54 S22-S17
      8 1279  918 (S05W12 L=204)   431    115        1 (1251/ 867-1309/ 950) -80.5 W10-W13 S08-S02
      9 1800  792 (S29W49 L=241)   867    173        0 (1762/ 727-1822/ 875)  53.1 W44-W51 S32-S23
     10 1512  666 (S27W21 L=214)   496    104        1 (1466/ 636-1549/ 713)   8.1 W17-W24 S28-S24
     11 1160  572 (S23E04 L=189)   847    280        1 (1028/ 505-1248/ 648)  10.6 E13-W03 S24-S19
     12 1443  419 (S43W12 L=204)  1335    370        0 (1297/ 326-1564/ 524)  13.5 E02-W23 S46-S38
     13  818  350 (S29E34 L=159)   407    106        0 ( 775/ 334- 877/ 368) -14.0 E36-E28 S29-S26
  

    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