Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1031 1776 (N52W34 L= 67)   395     60        1 (1007/1748-1053/1799)  59.0 W30-W36 N49-N54
      2 1181 1763 (N45W43 L= 76)   899     97        2 (1140/1727-1232/1790) -51.3 W38-W46 N42-N47
      3 1693 1535 (N14W64 L= 98)   898    119        1 (1640/1517-1752/1552) -11.3 W57-W72 N11-N16
      4 1126 1369 (N23W16 L= 49)   334     73        1 (1096/1341-1150/1401)  11.3 W13-W18 N21-N25
      5 1780 1313 (N00W59 L= 92)   638    131        1 (1751/1247-1801/1371)  74.1 W57-W61 S04-N03
      6  861 1162 (N19E06 L= 28)   679    203        1 ( 800/1093- 923/1244)  26.6 E10-W00 N16-N21
      7 1211 1144 (N08W14 L= 47)   264     51        0 (1187/1122-1228/1179)  45.0 W11-W15 N07-N10
      8  293  994 (N22E51 L=343)  1981    484        4 ( 257/ 759- 355/1151) -90.0 E54-E48 N07-N31
      9  478 1074 (N23E33 L=  1)   675    176        3 ( 392/1028- 548/1102) -28.6 E39-E28 N18-N25
     10 1872  937 (S24W58 L= 92)   633    103        0 (1841/ 899-1891/ 989)  26.6 W53-W61 S25-S21
     11 1555  859 (S17W27 L= 61)  1100    353        1 (1399/ 835-1722/ 882) -23.5 W17-W40 S22-S12
     12 1050  579 (S19E11 L= 23)   376    126        0 (1038/ 508-1065/ 638)  71.6 E12-E08 S24-S15
     13  705  523 (S14E33 L=  0)   839    256        0 ( 663/ 419- 748/ 628)  38.2 E40-E26 S19-S08
     14  610  434 (S18E44 L=350)   296     55        1 ( 584/ 421- 643/ 448) -45.0 E46-E41 S19-S16
     15  501  388 (S18E56 L=337)   375     31        1 ( 482/ 364- 518/ 414)  45.0 E58-E53 S20-S16
     16  584  333 (S24E53 L=340)   538    110        0 ( 539/ 312- 610/ 372)   7.1 E57-E48 S25-S21
     17  491  325 (S23E64 L=330)   998     95        5 ( 462/ 296- 532/ 362) -68.2 E71-E59 S25-S19
     18  783  314 (S30E39 L=354)   694    104        0 ( 749/ 270- 828/ 350)  11.3 E44-E34 S32-S27
  

    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