Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1408 1781 (N36W59 L=168)   502     79        1 (1378/1752-1448/1806) -14.0 W53-W63 N34-N37
      2 1034 1720 (N45W26 L=136)   251     22        0 (1016/1702-1053/1740) -71.6 W24-W27 N43-N47
      3  394 1611 (N60E43 L= 66)   595    108        1 ( 355/1563- 441/1641)  21.8 E48-E35 N56-N61
      4 1056 1384 (N24W12 L=121)   277     58        0 (1031/1362-1090/1401) -14.0 W09-W13 N22-N25
      5 1846 1338 (S01W66 L=176)  1120     80        2 (1798/1307-1893/1359) -12.1 W59-W73 S03-N01
      6  531 1128 (N22E27 L= 82)  1224    219        4 ( 423/1094- 645/1154) -17.4 E35-E19 N19-N26
      7 1080 1132 (N09W06 L=115)   175     84        0 (1033/1122-1128/1144) -21.8 W03-W08 N08-N10
      8 1479  861 (S16W22 L=132)   847    361        0 (1364/ 741-1575/ 989)  24.0 W12-W30 S19-S10
      9  790  907 (N04E16 L= 94)   403    159        0 ( 725/ 885- 845/ 943) -14.0 E19-E11 N02-N07
     10 1590  879 (S18W30 L=139)   382    122        1 (1543/ 832-1621/ 931)  29.7 W25-W32 S19-S15
     11  174  862 (N14E65 L= 45)  1505    229        2 (  94/ 833- 275/ 897)  -4.2 E79-E52 N12-N16
     12 1150  882 (S06W03 L=113)   191     96        0 (1104/ 866-1197/ 892) -36.9 W00-W05 S07-S04
     13  336  700 (N02E52 L= 57)  2360    556        3 ( 231/ 522- 415/ 891)  28.6 E76-E42 S06-N12
     14 1729  704 (S34W39 L=149)  1786    312        3 (1662/ 615-1824/ 849)  29.1 W32-W50 S38-S27
     15 1499  683 (S27W20 L=129)   250     68        0 (1478/ 650-1520/ 719)  36.9 W17-W21 S28-S25
     16  635  635 (S07E33 L= 77)   431    134        0 ( 594/ 585- 687/ 677)  20.6 E36-E28 S09-S05
     17 1043  272 (S42E24 L= 86)   418    115        0 ( 999/ 247-1084/ 319)  16.7 E28-E18 S43-S39
     18  892  186 (S44E42 L= 67)   947    154        2 ( 811/ 170- 962/ 211) -12.1 E49-E35 S46-S41
  

    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