Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1008 1883 (N58E37 L= 37)  1202    110        2 ( 954/1853-1056/1909)  50.2 E43-E32 N53-N62
      2  709 1557 (N24E34 L= 41)   511    125        1 ( 658/1534- 773/1587)   0.0 E37-E28 N22-N25
      3  592 1483 (N16E39 L= 36)   617    191        1 ( 565/1414- 609/1565) -59.0 E43-E35 N11-N21
      4 1753 1394 (N38W50 L=124)   403     63        0 (1734/1368-1774/1425) -36.9 W47-W52 N36-N39
      5 1329 1216 (N16W14 L= 89)  1185    441        4 (1182/1069-1438/1304) -21.8 W03-W23 N09-N19
      6 1526 1229 (N21W28 L=102)   234     43        1 (1504/1207-1547/1251)  71.6 W25-W29 N19-N22
      7 1554 1162 (N18W31 L=105)   464    102        2 (1508/1132-1600/1186)  33.7 W27-W33 N15-N19
      8  903  747 (S21E01 L= 73)   262    101        0 ( 859/ 723- 940/ 782)  59.0 E03-W00 S22-S17
      9 1631  594 (S14W52 L=126)   419     42        1 (1607/ 573-1658/ 616)   0.0 W49-W55 S15-S12
     10 1242  541 (S26W26 L=100)   890    249        1 (1156/ 511-1348/ 609)  42.0 W20-W30 S29-S19
     11  971  305 (S50W18 L= 93)   312     34        1 ( 948/ 282- 992/ 329) -26.6 W15-W21 S51-S47
     12 1247  237 (S47W48 L=122)  2516    261        5 (1146/ 207-1359/ 292)   6.8 W35-W60 S50-S42
  

    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