Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1740 1518 (N10W65 L=  4)   773    155        1 (1714/1471-1783/1554)  33.7 W61-W72 N07-N12
      2 1049 1432 (N27W13 L=312)   214     52        0 (1023/1409-1078/1451)   0.0 W11-W15 N25-N27
      3 1679 1405 (N07W52 L=350)   426    103        0 (1631/1384-1717/1425) -23.2 W49-W56 N05-N09
      4  179 1308 (N43E65 L=233)   702     74        0 ( 135/1288- 213/1323)   3.8 E74-E58 N41-N43
      5 1140 1164 (N10W10 L=308)  1346    527        2 ( 945/1039-1334/1247) -46.6 W00-W17 S02-N17
      6 1585 1186 (S01W37 L=335)   781    199        0 (1518/1132-1647/1226) -11.3 W30-W40 S04-N00
      7  616 1013 (N15E24 L=275)   954    199        3 ( 578/ 934- 668/1117)  85.2 E25-E21 N09-N21
      8 1613  782 (S25W29 L=327)  1948    406        0 (1515/ 620-1737/ 934)  26.6 W19-W41 S31-S19
      9  520  811 (N06E35 L=263)   644    133        2 ( 470/ 784- 588/ 850) -26.6 E38-E30 N03-N08
     10 1171  612 (S22E03 L=295)   651    157        1 (1122/ 546-1221/ 673)  32.0 E07-W01 S24-S18
     11  554  529 (S11E43 L=256)  1557    346        2 ( 482/ 470- 622/ 617) -66.0 E46-E37 S15-S04
     12 1435  500 (S37W11 L=309)   372    108        0 (1393/ 466-1473/ 543)  23.2 W07-W14 S38-S35
     13  410  468 (S11E58 L=241)   568     56        1 ( 386/ 435- 438/ 499) -76.0 E59-E55 S13-S09
     14 1021  453 (S27E17 L=281)   547     85        2 ( 984/ 425-1062/ 496)   0.0 E19-E13 S28-S24
     15  390  408 (S15E65 L=234)   505     85        1 ( 358/ 376- 412/ 438)  68.2 E68-E61 S17-S12
     16  551  362 (S21E53 L=246)  1255    174        2 ( 506/ 332- 603/ 411) -71.6 E57-E48 S23-S17
     17 1295  330 (S45E05 L=294)   311     54        0 (1276/ 302-1311/ 363)  45.0 E06-E02 S47-S42
     18  855  254 (S36E39 L=259)   712    104        2 ( 820/ 225- 892/ 289)   7.1 E42-E34 S38-S33
  

    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