Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1692 1594 (N26W59 L=  8)   455     80        0 (1667/1581-1715/1611)   0.0 W55-W61 N25-N27
      2 1758 1407 (N14W57 L=  7)   450     96        3 (1731/1386-1789/1424)   0.0 W54-W61 N12-N15
      3  185 1312 (N27E60 L=250)  3341    292        6 ( 136/1242- 225/1390)  45.0 E65-E56 N22-N32
      4  339 1282 (N23E44 L=266)   685    206        0 ( 257/1240- 402/1352)  22.6 E51-E38 N21-N26
      5 1732 1147 (S01W48 L=358)   499    107        0 (1710/1093-1760/1189)  78.7 W46-W51 S04-N01
      6 1210 1098 (N02W12 L=321)   231     39        0 (1191/1069-1226/1126) -90.0 W10-W12 N00-N03
      7  281  793 (S05E54 L=256)  1303    356        0 ( 240/ 664- 348/ 958) -86.6 E58-E50 S13-N04
      8  195  856 (S01E62 L=248)   431     85        3 ( 176/ 822- 215/ 891)  45.0 E63-E59 S02-N01
      9 1271  835 (S15W13 L=322)   192     38        0 (1257/ 801-1289/ 859)  71.6 W11-W13 S16-S13
     10  914  663 (S21E11 L=298)  1016    233        1 ( 828/ 566- 969/ 768)  42.5 E18-E06 S25-S14
     11  368  714 (S11E48 L=262)   470     99        0 ( 349/ 678- 391/ 757)  51.3 E50-E46 S13-S08
     12 1878  691 (S31W63 L= 12)   908     93        0 (1850/ 645-1902/ 730)  51.3 W60-W65 S33-S28
     13 1196  688 (S23W07 L=316)   206     64        3 (1162/ 651-1223/ 720)  45.0 W04-W08 S24-S20
     14  242  664 (S13E62 L=247)  1340    220        4 ( 197/ 616- 292/ 711) -10.3 E68-E56 S15-S09
     15 1447  629 (S30W24 L=333)   405     87        1 (1411/ 591-1487/ 663)  45.0 W21-W26 S31-S27
     16  304  600 (S17E58 L=252)  1262    237        4 ( 253/ 556- 335/ 649) -33.7 E62-E53 S20-S14
     17 1680  604 (S34W44 L=353)   524     91        0 (1664/ 558-1695/ 648) -71.6 W41-W45 S37-S31
     18  353  521 (S23E56 L=253)   800    178        2 ( 326/ 456- 372/ 587)  58.0 E59-E52 S26-S18
  

    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