Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  624 1738 (N31E48 L= 64)  1301    206        2 ( 541/1670- 696/1791) -16.9 E61-E37 N26-N34
      2  994 1677 (N33E17 L= 95)   769    115        1 ( 954/1622-1038/1725) -41.2 E20-E12 N29-N36
      3 1112 1486 (N22E04 L=108)   560    182        3 (1059/1441-1192/1551) -24.0 E08-W01 N19-N25
      4  934 1318 (N08E11 L=101)  1285    468        4 ( 901/1144- 980/1533) -61.4 E17-E04 S01-N21
      5 1663 1336 (N26W36 L=148)   794    121        1 (1631/1278-1698/1396) -54.5 W33-W39 N22-N29
      6 1289 1317 (N15W10 L=122)   651    227        3 (1233/1246-1379/1388) -26.6 W06-W16 N12-N19
      7 1271  836 (S13W18 L=130)   378    151        1 (1208/ 803-1345/ 869)  32.0 W14-W22 S15-S09
      8 1091  795 (S19W09 L=121)   265     38        0 (1067/ 776-1122/ 819)   0.0 W06-W10 S19-S17
      9  805  772 (S26E08 L=103)   426    115        0 ( 755/ 730- 857/ 809) -20.6 E12-E04 S27-S23
     10 1713  646 (S12W56 L=168)   389     52        2 (1693/ 618-1732/ 675) -36.9 W53-W58 S13-S10
     11  734  645 (S35E11 L=101)   332     73        0 ( 703/ 627- 775/ 671) -11.3 E13-E08 S36-S33
     12 1379  470 (S31W39 L=151)   346     87        1 (1334/ 455-1423/ 491)   0.0 W34-W42 S32-S29
     13 1540  439 (S28W53 L=165)   659    131        0 (1494/ 388-1575/ 490) -15.9 W46-W60 S30-S26
     14  840  449 (S46W01 L=113)   316     37        0 ( 815/ 428- 860/ 475)  71.6 E01-W02 S47-S44
     15 1326  231 (S47W56 L=168)   546     81        0 (1293/ 212-1368/ 258)  -4.1 W49-W63 S48-S46
     16  796  216 (S68W07 L=118)   920    171        1 ( 706/ 180- 868/ 253)  -2.0 E08-W20 S69-S65
  

    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