Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  512 1658 (N26E48 L=232)  1010    138        1 ( 468/1601- 568/1705)  82.9 E50-E44 N21-N29
      2 1105 1665 (N35E06 L=274)   436     62        1 (1078/1637-1141/1690)  45.0 E08-E03 N32-N37
      3 1191 1634 (N35W00 L=280)   393     76        0 (1167/1613-1228/1665) -26.6 E01-W03 N33-N36
      4 1092 1416 (N18E02 L=278)   599    110        2 (1066/1354-1125/1467) -81.9 E03-W00 N14-N21
      5 1345 1348 (N18W14 L=294)   772    276        1 (1286/1218-1381/1450)  86.4 W11-W16 N09-N25
      6  544 1386 (N08E36 L=244)   395     77        0 ( 515/1348- 565/1417) -31.0 E38-E33 N06-N10
      7 1396 1232 (N12W19 L=299)   522    146        0 (1350/1159-1443/1275)  66.0 W17-W21 N06-N15
      8 1698 1045 (N07W43 L=322)   391     54        1 (1680/1015-1720/1078) -45.0 W41-W44 N05-N08
      9 1798  987 (N06W52 L=332)   789     83        2 (1762/ 941-1831/1027)  50.2 W49-W55 N02-N08
     10  569  899 (S20E27 L=253)  4738   1329        6 ( 173/ 596-1002/1013)  13.7 E64-W00 S42-S08
     11 1296  668 (S22W23 L=303)   973    248        2 (1170/ 635-1416/ 706)  25.0 W15-W30 S26-S18
  

    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