Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1668 1602 (N24W68 L= 78)   751     83        3 (1631/1577-1694/1627) -45.0 W64-W73 N22-N26
      2  641 1446 (N40E19 L=352)  1161    221        0 ( 572/1399- 746/1508) -35.5 E24-E10 N35-N46
      3 1302 1341 (N19W24 L= 35)   257     96        0 (1288/1299-1322/1381)  78.7 W23-W26 N16-N21
      4  174 1289 (N36E72 L=299)  1584    131        4 ( 145/1241- 217/1341) -14.0 E85-E63 N32-N39
      5  315 1122 (N24E49 L=322)   507    116        0 ( 299/1086- 346/1168)  68.2 E50-E45 N22-N27
      6  226 1004 (N18E59 L=312)   423    122        0 ( 195/ 973- 261/1025)  15.9 E62-E54 N16-N19
      7  427  925 (N11E40 L=330)   389     96        0 ( 386/ 902- 473/ 941)   0.0 E43-E36 N10-N12
      8 1630  881 (S15W36 L= 46)   401    110        1 (1593/ 841-1658/ 932)  51.3 W32-W38 S16-S11
      9  222  812 (N06E62 L=308)   421     61        2 ( 197/ 789- 250/ 835) -21.8 E65-E59 N04-N07
     10  384  606 (S08E53 L=318)   942    142        1 ( 345/ 548- 427/ 665) -69.4 E54-E50 S11-S03
     11  598  603 (S11E35 L=335)   789    205        1 ( 507/ 541- 668/ 661)   8.7 E43-E30 S14-S07
     12  459  563 (S11E48 L=323)   334     63        1 ( 441/ 529- 487/ 590)  63.4 E49-E45 S13-S09
     13 1600  500 (S40W32 L= 43)   683     82        2 (1571/ 463-1630/ 534)  38.7 W29-W35 S43-S37
     14  671  458 (S21E36 L=335)   581    166        0 ( 600/ 435- 748/ 490) -29.1 E39-E30 S23-S18
     15 1136  396 (S35E07 L=  4)   252     53        0 (1106/ 383-1165/ 414)  11.3 E09-E04 S36-S34
     16  804  347 (S32E32 L=339)   414    139        0 ( 767/ 289- 836/ 406)  33.7 E38-E26 S35-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