Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1223 1637 (N34W32 L=136)   557    129        0 (1157/1612-1282/1656) -11.3 W26-W36 N31-N35
      2  853 1392 (N29E01 L=104)   508    125        1 ( 791/1365- 917/1416) -39.8 E04-W02 N26-N31
      3 1577 1324 (N06W41 L=145)   350    106        2 (1555/1269-1602/1363)  33.7 W38-W44 N03-N08
      4 1746 1234 (S03W51 L=155)  1390    305        2 (1653/1139-1848/1343) -56.3 W46-W58 S11-N04
      5  852 1254 (N21E04 L=100)   252     76        0 ( 829/1212- 886/1294)  31.0 E06-E01 N19-N22
      6 1024 1186 (N13W04 L=108)   179     59        0 ( 991/1167-1058/1209) -45.0 W02-W05 N11-N14
      7 1548 1047 (S08W31 L=135)   964    262        3 (1511/ 909-1570/1173)  71.6 W27-W33 S15-N00
      8  193 1025 (N23E57 L= 47)  2087    420        4 ( 101/ 981- 303/1091) -21.8 E67-E46 N18-N27
      9  783 1061 (N12E13 L= 91)   261    107        1 ( 721/1044- 834/1077)  -8.1 E16-E09 N10-N12
     10  397  955 (N15E40 L= 64)   426     98        1 ( 358/ 901- 425/ 993) -66.8 E42-E39 N10-N17
     11 1791  946 (S21W48 L=152)   341     88        0 (1764/ 927-1825/ 960)   0.0 W44-W51 S21-S19
     12 1214  801 (S13W05 L=109)   352    115        1 (1164/ 770-1265/ 848)  14.0 W01-W09 S14-S11
     13 1669  712 (S31W34 L=138)   286     51        1 (1643/ 687-1689/ 733)  45.0 W31-W36 S32-S30
     14 1387  578 (S31W11 L=115)   553     76        2 (1360/ 536-1411/ 613) -90.0 W09-W12 S34-S28
     15  351  568 (S07E57 L= 47)  1084    192        2 ( 285/ 537- 437/ 609) -26.6 E62-E50 S09-S03
     16  933  495 (S24E20 L= 84)   322    138        0 ( 903/ 439- 955/ 552)  45.0 E23-E16 S26-S20
     17 1108  463 (S31E10 L= 94)   229     66        1 (1084/ 426-1131/ 496)  26.6 E13-E07 S33-S29
  

    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