Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1320 1643 (N45W22 L=262)   493    101        0 (1278/1617-1368/1664)  12.5 W18-W27 N42-N46
      2 1590 1603 (N43W49 L=289)  1965    287        3 (1462/1564-1708/1636)   0.0 W34-W67 N39-N45
      3  556 1375 (N22E34 L=206)   272     55        1 ( 532/1356- 589/1396)  14.0 E35-E31 N20-N23
      4 1598 1131 (N11W38 L=278)  2291    425        5 (1519/ 950-1649/1298) -79.2 W32-W41 N00-N21
      5  833 1220 (N14E13 L=227)   423    193        0 ( 770/1150- 892/1295) -41.6 E17-E08 N09-N18
      6  274 1198 (N09E55 L=185)   896    101        0 ( 251/1157- 303/1254) -45.0 E57-E51 N06-N12
      7  700 1229 (N14E22 L=219)   252     71        0 ( 680/1187- 718/1254)  90.0 E23-E20 N10-N15
      8  230 1182 (N08E59 L=181)   360     87        0 ( 212/1134- 240/1215) -51.3 E62-E57 N04-N09
      9 1069 1126 (N09W02 L=242)   356    120        3 (1046/1078-1101/1183) -80.5 W00-W04 N06-N12
     10  896  997 (N00E08 L=232)   792    225        1 ( 845/ 905- 944/1081)  84.3 E10-E04 S05-N05
     11  480  863 (S10E35 L=205)   756    243        0 ( 437/ 772- 548/ 961) -57.5 E38-E30 S15-S04
     12 1124  695 (S17W08 L=248)  1843    609        3 (1024/ 469-1225/ 934)  73.8 W02-W12 S33-S02
     13  273  840 (S13E53 L=187)   328     86        0 ( 262/ 796- 285/ 884)  68.2 E54-E51 S15-S10
     14 1675  708 (S15W48 L=288)   817    253        1 (1650/ 613-1718/ 809) -61.7 W44-W54 S21-S08
     15 1756  709 (S15W56 L=296)  1129    139        2 (1720/ 649-1811/ 775) -63.4 W51-W62 S18-S10
     16 1377  287 (S47W39 L=279)   502     82        0 (1338/ 267-1409/ 321)  38.7 W35-W41 S49-S43
     17  811  242 (S55E17 L=223)  1922    186        6 ( 725/ 215- 911/ 281)   7.8 E28-E06 S59-S50
     18 1082  261 (S51W11 L=251)   499     84        1 (1036/ 243-1124/ 278)   0.0 W06-W14 S52-S49
  

    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