Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  963 1771 (N56E08 L=140)   494    102        0 ( 918/1752-1018/1786)   5.2 E12-E01 N53-N57
      2  501 1569 (N37E45 L=103)  2157    300        5 ( 377/1480- 646/1640) -12.2 E65-E28 N31-N42
      3 1665 1563 (N37W58 L=206)   591     56        0 (1638/1526-1691/1593)  26.6 W52-W63 N34-N39
      4  277 1482 (N30E67 L= 80)  1398    133        3 ( 248/1429- 320/1533) -12.8 E80-E58 N26-N33
      5 1721 1486 (N32W60 L=207)   469     49        1 (1701/1466-1744/1512) -71.6 W56-W62 N30-N33
      6  997 1439 (N29E02 L=145)   898    261        3 ( 879/1397-1115/1474) -11.3 E10-W05 N26-N31
      7 1766 1377 (N24W59 L=207)   465     44        0 (1748/1353-1785/1400)  31.0 W56-W62 N22-N25
      8  456 1360 (N23E41 L=106)   259     64        0 ( 420/1345- 492/1372)   9.5 E44-E38 N21-N23
      9  891 1114 (N09E08 L=139)   294     81        0 ( 860/1076- 915/1145) -45.0 E10-E06 N06-N10
     10  335  950 (S03E47 L=101)   382     80        1 ( 295/ 922- 372/ 982)  23.2 E50-E43 S04-N00
     11 1806  857 (S08W57 L=205)   727    153        0 (1791/ 778-1830/ 934) -59.0 W54-W61 S12-S02
     12  882  769 (S12E09 L=139)   452    139        3 ( 823/ 734- 943/ 804) -14.0 E12-E04 S14-S10
     13  741  623 (S22E19 L=129)   234     76        0 ( 706/ 611- 789/ 632)   9.5 E21-E15 S23-S21
     14 1464  552 (S27W32 L=180)   908    144        2 (1428/ 505-1517/ 623) -41.6 W28-W37 S30-S21
     15  538  414 (S38E41 L=107)  1101    150        1 ( 477/ 367- 588/ 488)  29.1 E50-E32 S42-S32
     16  844  207 (S58E20 L=128)   435     43        1 ( 823/ 184- 860/ 229)  45.0 E22-E16 S60-S54
  

    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