Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1106 1787 (N43E20 L=172)   969    208        1 (1060/1717-1158/1852) -32.7 E27-E13 N38-N48
      2  515 1690 (N21E54 L=137)   426     32        0 ( 491/1673- 541/1704)   0.0 E56-E52 N19-N21
      3  549 1604 (N15E47 L=145)   601    164        0 ( 532/1536- 565/1660) -41.2 E50-E42 N11-N19
      4 1247 1108 (N03W10 L=202)   234     72        0 (1209/1096-1292/1122)  11.3 W07-W12 N01-N04
      5 1773 1032 (N15W47 L=238)   332     84        0 (1759/ 984-1785/1078) -71.6 W45-W48 N11-N17
      6  244 1028 (S24E53 L=139)   731    231        1 ( 176/ 974- 315/1071)  30.3 E59-E46 S28-S20
      7 1422  906 (S03W26 L=217)  1258    288        1 (1379/ 771-1451/1068) -54.5 W19-W30 S10-N05
      8 1168  785 (S16W15 L=206)  1118    227        2 (1048/ 742-1273/ 838)  42.3 W08-W19 S21-S11
      9 1849  815 (N06W62 L=254)   485    167        0 (1814/ 802-1892/ 832)  10.3 W57-W68 N04-N07
     10  870  526 (S40W05 L=197)  1711    359        4 ( 703/ 495-1047/ 564)  15.6 E07-W18 S44-S35
     11 1454  461 (S26W45 L=237)  1309    274        0 (1317/ 435-1558/ 492)  31.0 W36-W52 S29-S20
  

    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