Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1101 1795 (N50W41 L= 36)   327     56        0 (1079/1773-1134/1813) -26.6 W37-W43 N48-N51
      2 1650 1549 (N16W60 L= 56)   616     61        2 (1623/1513-1674/1577)  68.2 W57-W63 N13-N18
      3 1065 1264 (N19W09 L=  4)  2328    528        3 ( 908/1115-1267/1430) -66.0 W03-W16 N05-N32
      4 1417 1365 (N14W33 L= 28)   289     99        1 (1392/1340-1444/1396) -63.4 W30-W35 N12-N16
      5 1479 1370 (N13W37 L= 33)   409     79        2 (1448/1332-1506/1394) -90.0 W35-W39 N10-N14
      6 1478 1287 (N08W34 L= 29)   481    167        2 (1439/1219-1534/1334)  24.0 W29-W38 N05-N10
      7  381 1197 (N33E39 L=316)   853    214        2 ( 295/1165- 472/1247) -25.2 E48-E31 N29-N37
      8  292  907 (N17E52 L=304)  1733    341        4 ( 171/ 861- 443/ 957) -17.1 E66-E39 N11-N21
      9 1166  799 (S09W02 L=357)   412    178        0 (1109/ 746-1241/ 842)   6.3 E02-W07 S10-S08
     10 1544  451 (S45W19 L= 15)   614     54        1 (1519/ 420-1569/ 478)  63.4 W17-W22 S48-S42
     11 1473  385 (S48W12 L=  7)   551     68        1 (1443/ 358-1508/ 414) -49.4 W09-W15 S51-S43
     12 1403  325 (S50W04 L=359)   639     84        1 (1378/ 299-1439/ 356)  18.4 W01-W07 S53-S46
  

    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