Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1455 1758 (N36W58 L=138)  1190     93        2 (1416/1715-1492/1788)  12.1 W51-W65 N33-N38
      2  786 1460 (N33E06 L= 75)  1733    379        1 ( 634/1383- 942/1518) -33.7 E16-W02 N24-N39
      3 1335 1428 (N19W28 L=108)  1108    341        1 (1182/1370-1466/1517)   2.5 W17-W40 N16-N21
      4  249 1389 (N41E57 L= 23)   795     69        1 ( 223/1361- 288/1424) -22.6 E62-E50 N38-N43
      5  261 1306 (N35E53 L= 27)  1038    101        3 ( 230/1259- 297/1353) -41.2 E57-E47 N31-N38
      6  239 1035 (N18E54 L= 26)  1399    215        5 ( 208/ 948- 287/1145) -77.0 E57-E50 N12-N25
      7 1093  799 (S12E01 L= 79)   398    137        0 (1074/ 730-1110/ 862)  60.3 E03-W01 S15-S08
      8 1596  766 (S25W31 L=112)   265     35        0 (1576/ 743-1616/ 788) -71.6 W29-W33 S26-S23
      9  827  614 (S17E21 L= 59)   325     78        0 ( 804/ 586- 859/ 648) -68.2 E21-E19 S19-S14
     10 1508  589 (S35W23 L=103)   285     50        0 (1480/ 572-1537/ 606) -11.3 W20-W25 S35-S33
     11  352  459 (S18E65 L= 15)  1226    136        3 ( 321/ 406- 392/ 516)  35.0 E71-E60 S21-S13
     12  498  410 (S23E53 L= 27)   900    168        2 ( 457/ 340- 535/ 495)  61.4 E58-E50 S28-S17
  

    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