Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1009 1837 (N56W16 L=114)   496     93        1 ( 966/1809-1048/1863)  -5.7 W11-W21 N53-N58
      2 1339 1664 (N36W34 L=132)  1350    362        1 (1171/1613-1486/1734)   5.4 W19-W51 N33-N39
      3 1571 1688 (N34W57 L=156)   611     63        0 (1537/1671-1609/1707) -15.9 W52-W61 N32-N35
      4  433 1455 (N34E38 L= 60)   326     36        1 ( 408/1430- 453/1481) -63.4 E40-E36 N31-N35
      5  277 1313 (N26E52 L= 47)   548     85        1 ( 235/1285- 315/1347) -26.6 E57-E47 N23-N29
      6  265 1065 (N11E51 L= 47)  1537    263        4 ( 229/ 947- 294/1173) -85.9 E55-E49 N03-N17
      7 1783 1045 (S08W52 L=151)  1087    256        1 (1733/ 928-1823/1154)  61.7 W46-W57 S14-S01
      8  203 1118 (N15E58 L= 41)   353     44        1 ( 184/1102- 229/1141) -11.3 E60-E55 N14-N16
      9  179 1018 (N09E61 L= 38)  1412    157        4 ( 133/ 969- 235/1076)  22.6 E67-E54 N06-N12
     10 1215  844 (S13W09 L=108)   282    120        0 (1175/ 806-1251/ 896)  45.0 W06-W12 S15-S10
     11 1601  730 (S25W36 L=134)  1453    365        0 (1537/ 572-1690/ 879)  62.2 W31-W42 S35-S16
     12 1867  810 (S23W63 L=162)  1088    157        0 (1818/ 761-1921/ 858) -26.6 W55-W73 S26-S19
     13 1161  624 (S26W04 L=102)   436    178        1 (1126/ 558-1188/ 717) -90.0 W02-W05 S30-S20
     14  588  620 (S19E34 L= 65)   427    138        1 ( 551/ 570- 642/ 659)  29.7 E37-E29 S22-S17
     15  694  620 (S21E27 L= 72)   265    115        0 ( 649/ 594- 742/ 654)  15.9 E29-E22 S21-S18
     16 1320  589 (S31W14 L=113)   527     99        0 (1269/ 557-1366/ 622)  23.2 W10-W17 S32-S29
  

    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