Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  730 1832 (N57E49 L=310)  1375    109        0 ( 689/1781- 782/1860) -66.0 E56-E43 N51-N61
      2  856 1616 (N40E19 L=340)   795    165        3 ( 788/1567- 929/1648) -13.0 E25-E12 N36-N41
      3  420 1504 (N27E52 L=307)   344     54        0 ( 386/1487- 446/1522) -12.5 E57-E48 N26-N28
      4 1340 1439 (N31W20 L= 18)   352    103        0 (1316/1400-1378/1481) -51.3 W17-W22 N28-N33
      5  241 1418 (N20E70 L=289)  1458    106        3 ( 210/1373- 273/1464) -14.0 E84-E63 N17-N22
      6 1120 1443 (N29W03 L=  2)   237     48        1 (1094/1429-1149/1462)  14.0 W01-W05 N28-N30
      7 1509 1392 (N29W33 L= 32)   370     68        0 (1478/1361-1538/1427) -45.0 W30-W34 N26-N30
      8  643 1329 (N18E28 L=331)   642    260        0 ( 549/1253- 732/1394) -23.6 E36-E20 N13-N21
      9 1726 1178 (N16W49 L= 48)   447     60        1 (1708/1138-1748/1207) -76.0 W47-W51 N13-N17
     10 1166 1082 (N07W08 L=  7)   470    127        0 (1133/1040-1213/1129)  45.0 W06-W11 N03-N09
     11 1712  648 (S17W55 L= 54)   391     34        2 (1692/ 627-1732/ 671)  56.3 W52-W57 S18-S15
     12 1481  556 (S25W37 L= 36)   364     89        1 (1440/ 539-1528/ 580)  -7.1 W32-W40 S25-S23
  

    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