Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1197 1791 (N49W01 L= 10)   598    102        2 (1164/1755-1226/1826) -63.4 E02-W03 N45-N52
      2 1059 1806 (N47E12 L=358)   303     39        0 (1037/1784-1084/1823)  63.4 E13-E09 N45-N49
      3 1502 1590 (N38W27 L= 36)   504     85        0 (1476/1554-1537/1614)  45.0 W24-W30 N34-N39
      4  498 1437 (N13E40 L=330)   250     54        1 ( 468/1417- 522/1457) -11.3 E42-E37 N12-N14
      5 1883 1196 (N19W60 L= 69)   955    213        6 (1847/1164-1921/1229) -18.4 W55-W65 N16-N20
      6 1383 1157 (N07W20 L= 29)   215     49        2 (1367/1123-1400/1185) -63.4 W18-W20 N04-N08
      7  350 1138 (S06E45 L=324)   397     62        0 ( 330/1108- 370/1170) -45.0 E47-E43 S07-S04
      8 1888 1133 (N15W61 L= 70)   475    131        4 (1863/1105-1915/1163) -15.9 W57-W64 N13-N16
      9  744 1097 (S05E18 L=352)   365    141        1 ( 685/1032- 796/1150) -36.9 E21-E13 S08-S02
     10  207 1044 (S12E58 L=311)   433    137        3 ( 176/1016- 246/1076)   7.1 E61-E53 S14-S10
     11  709  942 (S15E18 L=351)   547    155        0 ( 662/ 887- 789/ 984)  45.0 E20-E13 S18-S11
     12  142  891 (S22E65 L=305)  1327    210        4 ( 110/ 850- 185/ 944)  29.1 E69-E59 S23-S18
     13 1333  878 (S10W21 L= 30)   379    106        1 (1288/ 830-1366/ 932)  74.1 W18-W21 S13-S06
     14  309  891 (S21E47 L=322)   277     80        1 ( 296/ 856- 324/ 930) -78.7 E48-E46 S23-S18
     15 1019  819 (S18W02 L= 12)   188     68        0 ( 978/ 809-1058/ 830)  14.0 W00-W04 S18-S17
     16 1658  807 (S09W44 L= 54)   246     31        0 (1634/ 793-1679/ 822)  18.4 W42-W46 S09-S07
     17  492  701 (S31E33 L=337)   480     98        1 ( 460/ 668- 537/ 741) -26.6 E35-E29 S33-S29
     18 1797  563 (S20W66 L= 76)   594     72        3 (1778/ 535-1818/ 584)  76.0 W64-W69 S21-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