Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1542 1691 (N30W59 L= 72)   563     84        0 (1522/1669-1564/1711)  11.3 W56-W61 N28-N31
      2 1824 1033 (S14W54 L= 68)  1443    161        3 (1773/ 977-1858/1097)  16.7 W48-W58 S17-S10
      3  240 1047 (N24E55 L=319)   875    166        1 ( 191/1001- 287/1085) -36.9 E59-E50 N20-N26
      4 1169  825 (S07W04 L= 17)   336    108        0 (1115/ 793-1224/ 846)  -9.5 W00-W06 S08-S06
      5  465  712 (N00E43 L=331)  1044    200        3 ( 437/ 618- 492/ 803)  68.2 E45-E38 S05-N05
      6 1345  652 (S22W10 L= 24)   624    231        0 (1299/ 543-1366/ 762)  58.4 W05-W13 S28-S15
      7  988  699 (S11E10 L=  4)   192     64        1 ( 965/ 671-1004/ 740)  63.4 E11-E08 S12-S08
      8 1457  596 (S30W17 L= 31)   648    145        1 (1411/ 535-1508/ 663)  36.9 W12-W21 S32-S26
      9  963  422 (S28E20 L=354)   608    216        0 ( 883/ 362-1057/ 491)  17.4 E27-E11 S30-S25
  

    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